Low Mass Dimuons Produced in Relativistic Nuclear Collisions 
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The NA60 experiment has measured low-mass muon pair production in In-In coUisions at 158 A 
GeV with unprecedented precision. We show that this data is reproduced very well by a dynamical 
model with parameters scaled from fits to measurements of hadronic transverse mass spectra and 
Hanbury-Brown and Twiss correlations in Pb-Pb and Pb-Au collisions at the same energy. The 
data is consistent with in-medium properties of p and ci;-mesons at finite temperature and density 
as deduced from empirical forward-scattering amplitudes. Inclusion of the vacuum decay of the 
p-meson after freeze-out is necessary for an understanding of the mass and transverse momentum 
spectrum of dimuons with M < 0.9 GeV/c^. 

PACS numbers: 25.75.-q,25.75.Gz 



I The main goal of the relativistic nuclear collision pro- 
gram is to produce and study strongly interacting matter 

] at high temperature and density. It is hoped that exotic 
many-body effects may be uncovered, one of them being a 

] quark-gluon plasma (QGP), a state where hadronic mat- 
ter exhibits partonic behavior. However, other interest- 

I ing phases may also manifest themselves. In this context, 
electromagnetic observables - real and virtual photons - 

I constitute a privileged class of probes because of the near 
absence of final state effects. The radiation will travel 

I essentially unscathed from its production point to the 
detectors. As the system expands and cools, a quantita- 

" tive understanding of the net electromagnetic spectrum 
requires a detailed understanding of the local emissivity 

] as well as knowledge of the space-time evolution of the 
radiating matter. 

] The recent NA60 experiment at the CERN Super- 
Proton Synchrotron (SPS) has measured the production 

I of low-mass muon pairs in In-In collisions at 158 A GeV. 

' In this experiment the spectra of invariant mass (M) in 
the region M < 1.5 GeV/c^, and of transverse momen- 
tum (pt), were obtained with unprecedented precision 
[l[. While the invariant mass spectrum is essential to 
characterize in-medium changes to the electromagnetic 
current-current correlation function, the transverse mo- 
mentum spectrum of dilcptons is especially sensitive to 
the interplay between production processes and collective 
transverse flow. A simultaneous description of pr and M 
spectra therefore constitutes a stringent test of the dy- 
namical evolution of the produced matter and of our un- 
derstanding of in-medium modifications of vector mesons 
as revealed through thermal dilepton production. In this 
paper we pursue a theoretical interpretation of the re- 
cent NA60 data that involve folding microscopic dilepton 
emission rates with a dynamical evolution model. The 



local pair production rate may be written as [4] 
^ = ^P(M)W,-)/b(.o,T) 

where a is the electromagnetic fine-structure constant, 
T is temperature, and M, go and q are the dimuon 
mass, energy and momentum, respectively. The func- 
tion P{M) accounts for the phase-space reduction due 
to the finite rest mass of the muon and does not depend 
on the medium's properties. The function /b is the Bose- 
Einstein distribution. The virtual photon spectral func- 
tion (averaged over polarizations), R{M, q), is directly re- 
lated to the retarded in-medium electromagnetic current- 
current correlator IIJ^^^ via R = — {Air / M'^)ItiiIV'^^ ^. 

In the low invariant mass region, M < \ GeV/c^, 
the in-medium modifications of this correlator in the 
hadronic phase are directly related to spectral proper- 
ties of light vector mesons {p, oj and (j)) through vec- 
tor meson dominance (VMD) Of these, the con- 
tribution from the p-mcson is the largest. Our focus 
here is on contributions from the p and lo mesons; the 
modifications of the (/)-meson will be studied elsewhere. 
The in-medium vector meson spectral densities are eval- 
uated in an approach where it is assumed that the vector- 
isovector and vector-isoscalar fields are modified mainly 
through scattering from nuclcons and pions in the heat 
bath [J] . One can infer the finite temperature and density 
dependence of the spectral functions in the single scat- 
tering limit by extracting the leading term in the self- 
energy expansion. More specifically, the self-energy is 
related to forward scattering amplitudes which are eval- 
uated in a two-component model that involves the exci- 
tation of s-channel resonances upon a background - dual 
to the Pomeron - which becomes important especially at 
high energies 0, 0] ■ Alternatively, effective hadronic La- 
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FIG. 1: The imaginary part of the p-meson propagator for 
temperature T = 150 MeV and baryonic density 0.5 and 1 
(in units of nuclear matter density). The imaginary part of 
the vacuum propagator is also shown. 



grangian 0, 0] or chiral reduction techniques Q may be 
used. The results from both methods yield a negative 
contribution to the pole mass from interactions with pi- 
ous and a positive contribution from interactions with 
nucleons; the net deviation from the vacuum mass be- 
ing small (at most a few tens of MeV) at temperatures 
and densities such as those probed in 158 A GeV In- 
In collisions. The spectral width, however, is increased 
considerably owing to interactions with the medium [^; 
these features are shown in Figure [T] Collisions with nu- 
cleons are the dominant effect but collisions with pious 
also contribute 0, Q . 

Two additional sources are emission from a thermal- 
ized partonic phase and from four-pion annihilation pro- 
cesses. A deconfined phase existing at early times and 
high temperatures is modeled by means of a quasipar- 
ticle model @. Previous studies have made clear that 
quantifying dimuon emission in the invariant mass re- 
gion M > 1 GcV / (? requires the inclusion of four-pion 
states [lO[. In the present work the required matrix el- 
ements were evaluated in the Lagrangian approach of 
Ref- [lit where the inverse process, e'^e~ annihilation 
in the vacuum into four charged pious, was studied. The 
model was shown to provide reasonable agreement with 
the measured cross sections. In addition, we included the 
annihilation of two neutral and two charged pious within 
the same framework supplemented by intermediate states 
containing the uj and hi mesons. 

Two contributions other than radiation from the ther- 
malized expanding medium are essential to the under- 
standing of the experimentally observed spectra in the 
mass region M < 1.5 GeV/c^: the dileptons emitted 
from hadronic decays after the system has blown apart, 
and those from correlated charm decays {D, D). The 



late hadronic decay contribution (for example, the spec- 
tral profile contributions from vacuum (p and uj decay) 
has been subtracted in the M and pt distribution by 
NA60, but not that of the freeze-out p's. The spectra 
analyzed by NA60 therefore reflect not only the excess 
dimuon spectrum, but also include a substantial contri- 
bution from the decay of vacuum p-mesons outside of the 
medium after hadronic freeze-out. Since those p-mesons 
are strongly influenced by transverse flow, it is essential 
to go beyond the statistical model formulation and cal- 
culate their momentum spectrum after decoupling using 
the Cooper-Frye formula [11] • This again stresses the 
need for a realistic description of the spatial and tempo- 
ral evolution of the system. 

A detailed description of the fireball evolution model 
for nucleus- nucleus collisions at the CERN SPS can be 
found in [l^ . ITsI . The main assumption is that an equili- 
brated system is formed a short time tq after the nuclear 
impact. The fireball subsequently expands isentropically 
until mean free paths exceed the system size and particles 
free-stream to the detectors. This is assumed to occur at 
a freeze-out time (or temperature) that is the same for 
all hadronic species. The entropy density s(r, rjs, r) is de- 
scribed by the product of two Woods-Saxon distributions 
s = NR(r,T)H{ris,T) that depend on the spacetime ra- 
pidity rjs = ill! (^f^t^ s-iid the transverse-plane radius r. 
The is a normalization constant. The Woods-Saxon 
profiles R{r,T) = (1 -I- exp [(r — i?c(T)) /dwsD^^ and 
H{r]s,T) = (14- exp [(|r/s| - i7c(T)) /?7ws])~^ are charac- 
terized by the thickness parameters d^s and yy^s and by 
the size of the emitting zone as a function of proper time 
via Rc{t) and Hc{t). The latter two functions are cal- 
culated under the assumption of constant radial and lon- 
gitudinal acceleration, respectively. This translates into 
the rapidity of the fireball front being ?7s(t) = ?]o + driT. 
The parameter 2r/o is the initial size occupied by the fire- 
ball and a^i is a longitudinal expansion parameter. The 
initial radial extension Rc{Q) is determined in a calcula- 
tion of the initial density profile using the Glauber model. 
Transverse flow is described best if transverse rapidity pT 
scales like -Jr [iJl- The accelerated longitudinal expan- 
sion driven by strong initial compression implies that in 
general space-time and momentum rapidity are not the 
same; their mismatch A77 can be seen as a characteri- 
zation of how much the solution departs from the ideal 
Bjorken [ist scenario. The parameters rg, a^, flj;, ^yo, to- 
gether with the decoupling temperature Tf, set the scale 
of the spacetime evolution, and d^s and ?7ws specify the 
details of the entropy density. 

The parameters listed above should in principle be fit 
to hadronic data and then used to predict the dimuon 
spectra. However, the necessary hadronic data for In- 
In collisions is not yet available; fortunately, the data for 
Pb-Pb and Pb-Au collisions at the same beam energy are 
0, III- The same theoretical framework as described 
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here has proven successful in the description of pho- 
ton and dilepton emission and charmonium supression 
in those collisions _13;]. For the collisions of In on In at 
the SPS, the total entropy in semi-central In-In collisions 
at 77 = 3.8 measured by NA60 [l[ is obtained from that 
in peripheral (30%) Pb-Au coUisions with 2.1 < ?; < 2.55 
measured by CERES [3 by multiplying by the ratio 
of charged particle rapidity densities dNch/drj measured 
in both experiments. The number of participant baryons 
and initial spatial extent are obtained via geometrical nu- 
clear overlap calculations, while rjo is determined under 
the assumption that stopping power scales approximately 
with the number of binary collisions per participant. The 
electromagnetic emission near midrapidity turns out to 
be quite insensitive to changes in the values of 77^3 and 
dvis- these are therefore not modified. The parameters of 
the accelerated expansion a± and well as the equi- 

libration time are assumed to be primarily determined 
by the incident system energy so they are kept as in Pb- 
Pb collisions. The largest uncertainty is the choice of 
the decoupling temperature. Because the In-In system 
is smaller than Pb-Pb, a higher decoupling temperature 
is expected. However, to give an unambiguous answer 
would only be possible with simultaneous measurements 
of HBT correlations and transverse mass spectra. Here 
we choose Tf = 130 McV. As alluded to earlier, kinetic 
equilibrium of all processes is assumed until this univer- 
sal decoupling temperature is reached. The fact that the 
four-pion processes contribute all the way down to Tf 
translates into an upper limit to their contribution, as 
microscopic descriptions of the dynamics (see, for exam- 
ple, [1^) suggest a sequential decoupling of the different 
channels. The equation of state in the hadronic phase 
and the chemical potentials (It^ , fiK are inferred from sta- 
tistical model calculations as described in HQ. The 
resulting fireball has a peak temperature of about 250 
MeV and a lifetime of about 7.5 fm/c. 

The thermal contributions have to be calculated by 
folding the space-time evolution of the expanding mat- 
ter with the thermal rates. For the mass spectra this 
amounts to 
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where A represents the detector acceptance of NA60. 
The space-time evolution enters via the dependence of 
the thermal rates on temperature, baryon and pion chem- 
ical potential, and energy and momentum of the decaying 
virtual photon in the local rest frame. The equilibrium 
thermal rates have been augmented by appropiate fugac- 
ity factors exp(n/U7r/r) with n = 2,3,4 for the thermal 
p, Lu, and four-pion annihilation contributions. 

A comparison of a theoretical calculation of the in- 
variant mass spectra, with its different components, with 
NA60 measurements in semi-central In-In collisions at 
the SPS is presented in Fig. ^ Consider the all pr- 
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FIG. 2: Mass spectrum in semi-central In-In collisions at SPS 
[l| compared to theory. The lower panel shows an integration 
over all pr, the upper panel one over high transverse momenta 
1 GeV /c < pt < 2 GeV / c. Partial contributions arise from p 
decays in vacuum after freeze-out, thermal in-medium p and 
u decays, radiation from a thermalized QGP and from ther- 
mal four-pion annihiliation, and from correlated open charm 
decay. 



data in the lower panel first. The region below 1 GeV 
and at masses smaller than the vacuum mass of the p 
meson clearly demonstrates a considerable in-mcdium 
broadening of the p and lo mesons. The decays of the 
vacuum p after freeze-out are important, and in certain 
mass ranges are of the same order as the in-mcdium p 
and LO meson decay contributions. The contribution from 
the QGP phase becomes more important and eventually 
dominant with higher M . We find that four-pion annihi- 
lation processes are subdominant even for masses above 
M > 1.25 GeV/c2. This differs from the findings in [2lj 
where the four-pion annihilation was derived in the soft 
pion limit, assuming chiral mixing. 

Even though the spectrum is integrated over all pr, 
and mass is Lorentz invariant, a successful theoretical 
understanding of these spectra still requires a detailed 
understanding of the px spectra since the acceptance A 
is transverse momentum-dependent. Additional informa- 
tion is obtained if one compares the theoretical predic- 
tion for momentum cuts of the mass spectrum with the 
experimental data. Since the acceptance restricts con- 
siderably the information that can be obtained from the 
mass spectrum at low transverse momenta, this contri- 
bution is not shown even though our model does provide 
a good description of the spectrum in this window of 

< Pt < 0.5 GeV/c. At higher transverse momenta, 

1 GeV/c < Pt < 2 GeV/c, the relative contribution of 
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production as measured in 158 A GeV In-In collisions at 
the CERN SPS reflects substantial in-medium broaden- 
ing of the p meson spectral function in the hot and dense 
nuclear medium. Furthermore, we found that at higher 
invariant masses thermal radiation with T > 170 MeV 
dominates over four pion annihiliation processes. This 
is especially relevant to a theoretical understanding of 
intermediate mass dimuon production. 
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FIG. 3: Transverse momentum spectra as obtained from the- 
ory for semi-central In-In collision. The acceptance-corrected 
data are from NA60 Q]. 



dimuons from vacuum p-dccays is relatively enhanced as 
the emission then occurs at a later stage of the evolution 
where considerable transverse flow has already built up. 
The robustness of our theoretical understanding of the 
collision can be further assessed from a comparison with 
the acceptance corrected pt spectra for different mass 
windows: see also IJ]. Figure [3] shows the pt spectra in 



three different mass windows for semi-central collisions. 
(The data are averaged over different centrality classes 
excluding peripheral collisions We also performed 



this averaging in [IJ, [23] and found that differences be- 
tween the averaged data and the semi-central collision 
data are small.) 

The mass region 0.4 GeV/c^ < M < 0.6 GeV/c^ and 
the p-like mass region 0.6 GeV/c^ < Af < 0.9 GeV/c^ 
receive most of their contribution from the late hadronic 
stages, namely, from decays of in-medium vector mesons 
and the vacuum p mesons after freeze-out. In those stages 
considerable flow has already been built up which im- 
plies that the blueshift of the spectra by flow is large. 
The difference between the p-like and the lower mass re- 
gion in the spectra is caused by the different contribu- 
tions of the in-medium vector mesons and the vacuum 
p. The latter receives the maximum ffow and predomi- 
nantly contributes to the p-like region for momenta above 
~ 1 GeV/c. The slope of the transverse momentum spec- 
trum in the mass region 1.0 GeV/c^ < M < 1.4 GeV/c^ 
is dominated by contributions from the early QGP phase 
where flow has not yet built up while a contribution 
from four- pion annihilation processes is subdominant. If 
four-pion annihilation processes were more substantial 
this would result in considerable hardening of the pT- 
spectrum which is not observed 0,0, 22|. 

In conclusion, wc have shown that low mass dimuon 
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